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Abstract 

(N 

,— H We present the first of two articles on the small volume fraction limit of a nonlocal 

^ ^ Cahn-Hilliard functional introduced to model microphase separation of diblock copoly- 

mers. Here we focus attention on the sharp-interface version of the functional and con- 
sider a limit in which the volume fraction tends to zero but the number of minority phases 
(called particles) remains 0(1). Using the language of F-convergence, we focus on two 
I— I levels of this convergence, and derive first and second order effective energies, whose en- 

ergy landscapes are simpler and more transparent. These limiting energies are only finite 
on weighted sums of delta functions, corresponding to the concentration of mass into 
(-H 'point particles'. At the highest level, the effective energy is entirely local and contains 

information about the structure of each particle but no information about their spatial 
distribution. At the next level we encounter a Coulomb- like interaction between the par- 
ticles, which is responsible for the pattern formation. We present the results here in both 
three and two dimensions. 
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1 Introduction 

This paper and its companion paper [11] are concerned with asymptotic properties of two 
energy functionals. In either case, the order parameter u is defined on the fiat torus T" = 
R'^/Z", i.e. the square [—5, 5]"" with periodic boundary conditions, and has two preferred 
states M = and u = 1. We will be concerned with both n = 2 and n = 3. The nonlocal 
Cahn-Hilliard functional is defined on /f^(R") and is given by 

S^iu) ■■= e I \\7u\'^dx + - I u^{l-u^)dx + cr ||n - f nH^^i.^^-,. (1.1) 

Its sharp interface limit (in the sense of L-convergence), defined on 5y(T"; {0, 1}) (charac- 
teristic functions of finite perimeter) , is given by [25] 



£{u) := / |Vu| + 7 ll-u - f -ull^-i/rpnA. (1.2) 



In both cases we wish to explore the behavior of these functionals, including the structure of 
their minimizers, in the limit of small volume fraction Jrj.„ u. The present article addresses the 



sharp interface functional ( 1.2 ); the diffuse- interface functional S"^ is treated in the companion 
article ITTI. 



1.1 The diblock copolymer problem 

The minimization of these nonlocal perturbations of standard perimeter problems are natural 
model problems for pattern formation induced by competing short and long-range inter- 
actions [33]. However, these energies have been introduced to the mathematics literature 
because of their connection to a model for microphase separation of diblock copolymers 0. 

A diblock copolymer is a linear-chain molecule consisting of two sub-chains joined cova- 
lently to each other. One of the sub-chains is made of monomers of type A and the 
other consists of Nb monomers of type B. Below a critical temperature, even a weak repul- 
sion between unlike monomers A and B induces a strong repulsion between the sub-chains, 
causing the sub-chains to segregate. A macroscopic segregation where the sub-chains detach 
from one another cannot occur because the chains are chemically bonded. Rather, a phase 
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separation on a mesoscopic scale with A and B-rich domains emerges. Depending on the ma- 
terial properties of the diblock macromolecules, the observed mesoscopic domains are highly 
regular periodic structures including lamellae, spheres, cylindrical tubes, and double-gyroids 
(see for example [6l). 

A density- functional theory, first proposed by Ohta and Kawasaki |2T], gives rise to a 
nonlocal free energy [20' in which the Cahn-Hilliard free energy is augmented by a long-range 
interaction term, which is associated with the connectivity of the sub-chains in the diblock 
copolymer macromoleculej^ 
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iViipdx + / u'^{l-u'^)dx + ^ ||u-M||^_l,T„^. (1.3) 
Often this energy is minimized under a mass or volume constraint 

u = M. (1.4) 



/ 



Here u represents the relative monomer density, with u = corresponding to a pure-j4 region 
and u = 1 to a pure-i? region; the interpretation of M is therefore the relative abundance 
of the A-parts of the molecules, or equivalently the volume fraction of the A-region. The 
constraint of fixed average M reflects that in an experiment the composition of the molecules 



is part of the preparation and does not change during the course of the experiment. In (1.3) 
the incentive for pattern formation is clear: the first term penalizes oscillation, the second 
term favors separation into regions oi u = and u = 1, and the third favors rapid oscillation. 



Under the mass constraint (1.4) the three can not vanish simultaneously, and the net effect 



is to set a fine scale structure depending on e, cr and M. Functional (1.1 1 is simply a rescaled 



version of (1.3) with the choice of o" = £7. Its sharp-interface (strong-segregation) limit, in 



the sense of F-convergence, is then given by (1.2) ^ 



1.2 Small volume fraction regime of the diblock copolymer problem 

The precise geometry of the phase distributions (i.e. the information contained in a minimizer 



of (1.3)) depends largely on the volume fraction M. In fact, as explained in [lOj, the two 
natural parameters controlling the phase diagram are £^/o^ and M. When the combination 
£y/a is small and M is close to or 1, numerical experiments [TU] and experimental obser- 
vations [6] reveal structures resembling small well- separated spherical regions of the minority 
phase. We often refer to such small regions as particles, and they are the central objects of 
study of this paper. 

Since we are interested in a regime of small volume fraction, it seems natural to seek 
asymptotic results. It is the purpose of this article and its companion article [H] to give a 
rigorous asymptotic description of the energy in a limit wherein the volume fraction tends to 
zero but where the number of particles in a minimizer remains 0(1). That is, we examine 
the limit where minimizers converge to weighted Dirac delta point measures and seek effec- 
tive energetic descriptions for their positioning and local structure. Physically, our regime 
corresponds to diblock copolymers of very small molecular weight (ratio of B monomers to 



^See |T2] for a derivation and the relationship to the physical material parameters and basic models for 
inhomogeneous polymers. Usually the wells are taken to be ±1 representing pure phases of A and _B-rich 
regions. For convenience, we have rescaled to wells at and 1. 
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Figure 1: Top: an AB diblock copolymer macromolecule of minority A composition. Bottom: 
2D schematic of two possible physical scenarios for the regime considered in this article. Left: 
microphase separation of very long diblock copolymers with minority A composition. Right: 
phase separation in a mixture/blend of diblock copolymers and homopolymers of another 
monomer species having relatively weak interactions with the A and B monomers. 



A), and we envisage either a melt of such diblock copolymers (c/. Figure [T] bottom left) or a 
mixture/blencj^of diblocks with homopolymers of type A {cf. Figure [l] bottom right). 

This regime is captured by the introduction of a small parameter r/ and the appropriate 
rescaling of the free energy. To this end, we fix a mass parameter M reflecting the total amount 
of minority phase mass in the limit of delta measures. We introduce a small coefficient to M, 
and consider phase distributions u such that 

u = r/"M, (1.5) 

where n is either 2 or 3. We rescale u as follows: 



u 



V := — , (1.6) 
so that the new preferred values of v are and l/ry". We now write our free energy (either 



(1.1) or (1.2)) in terms of v and rescale in rj so that the minimum of the free energy remains 



0(1) as r/ ^ 0. In this article, we focus our attention on the sharp interface functional (1.2): 



that is, we assume that we have already passed to the limit as e — > 0, and therefore consider 



the small- volume-fraction asymptotics of ( 1.2 ). In [TTl we will show how to extend the results 



of this paper to the diffuse-interface functional (1.1), via a diagonal argument with a suitable 
slaving of e to r]. 

In Section [s] we consider a collection of small particles, determine the scaling of the H~^- 
norm, and choose an appropriate scaling of 7 in terms of rj so as to capture a nontrivial limit 
as r] tends to 0. This analysis yields 



'r?E2d(r;) if n = 2 



r?2E3d(v) ifn = 3. 



A similar nonlocal Cahn-Hilliard-like functional models a blend of diblocks and homopolymers [13j . 
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where 



E^^^^) := V j J'^v\ + \^ogr]\~^\\v-f^,v\\l_,^^,^ defined for v £ BV {T^; {0,1/r]^}) 

(1.7) 



and 



E^^\v) ■= V j\ |Vw| +T]\\v- ^^||^-i(T3) defined for v e BV (T^; {0, l/rf^]) . (1.8) 

In both cases, E^^(u), E^{v) remain 0(1) as ry — > 0. 

The aim of this paper is to describe the behavior of these two energies in the limit ry — > 0. 
This will be done in terms of a T -asymptotic expansion for ^^{v) and E^'^(t>). That is, we 
characterize the first and second term in the expansion of, for example, E^^ of the form 

E^*^ = Eq*^ + 7/ Fg'^ + higher order terms. 

Our main results characterize these first- and second-order functionals Eq'^ , Fq'^ (respec- 
tively Eq*^, Fq*^) and show that: 

• At the highest level, the effective energy is entirely local, i.e., the energy focuses sep- 
arately on the energy of each particle, and is blind to the spatial distribution of the 
particles. The effective energy contains information about the local structure of the 



small particles. This is presented in three and two dimensions by Theorems 4.2 and 6.1 
respectively. 

At the next level, we see a Coulomb-like interaction between the particles. It is this 
latter part of the energy which we expect enforces a periodic array of particles]^ This 



is presented in three and two dimensions by Theorems 4.4 and 6.4 respectively. 



The paper is organized as follows. Section [2] contains some basic definitions. In Section [3] 
we introduce the small parameter ry, and begin with an analysis of the small-ry behavior of the 
norm via the basic properties of the fundamental solution of the Laplacian in three and 
two dimensions. We then determine the correct rescalings in dimensions two and three, and 
arrive at ([LT]) and In Section |4] we state the T-convergence results in three dimensions, 

together with some properties of the T-limits. The proofs of the three-dimensional results are 
given in Section [5] In Section [6] we state the analogous results in two dimensions and describe 
the modifications in the proofs. We conclude the paper with a discussion of our results in 
Section [7l 



2 Some definitions and notation 

Throughout this article, we use T" = R"/Z" to denote the n-dimensional flat torus of unit 
volume. For the use of convolution we note that T" is an additive group, with neutral element 
G T" (the 'origin' of T"). For v G BV{T''] {0, 1}) we denote by 

''Proving this is a non-trivial matter; see Section 7 
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the total variation measure evaluated on T", i.e. HVnlKT") Since v is the characteristic 
function of some set A, it is simply a notion of its perimeter. Let X denote the space of 
Radon measures on T". For /ir;,/^ € X , fi^j fj, denotes weak-* measure convergence, i.e. 



for all / G C(T"). We use the same notation for functions, i.e. when writing t;^ ^ vq, we 
interpret and vq as measures whenever necessary. 

We introduce the Green's function Gj^n for —A in dimension n on T*^. It is the solution 

of 



-AGt" =5-1, with / Gt" = 0, 

where 6 is the Dirac delta function at the origin. In two dimensions, the Green's function 
G^a satisfies 

Gt,2{x) = -^log\x\ + g^^\x) (2.1) 

for all X = (xi, X2) G with max{|xi|, \x2\} < 1/2, where the function g^'^^ is continuous on 
[-l/2,l/2]2 and C°° in a neighborhood of the origin. In three dimensions, we have 

Gt^{x) = ^ + (2.2) 

for all X = {xi,X2,X3) G with max{|xi|, \x2\, ja^sl} < 1/2, where the function g^^^ is again 
continuous on [-1/2, 1/2] 3 and smooth in a neighbourhood of the origin. 
For ji e X such that /^(T") = 0, we may solve 

—Av = 

in the sense of distributions on T'^. If ?; G H^{T^), then fi G i?-^(T"), and 

ll/^ll//~i(T") / I^^P dx. 

In particular, if u G L2(T") then (n - fn) G F"^(T") and 

11^ ~ f ^ll?/-i(T") = / / u{x)u{y) GTn{x - y) dx dy. 

Note that on the right-hand side we may write the function u rather than its zero-average 
version u — f u, since the function Gt" itself is chosen to have zero average. 

We will also need an expression for the H^^ norm of the characteristic function of a set 
of finite perimeter on all of R'^. To this end, let / be such a function and define 

Ih-i(R^) = I'^^l'^dx, 
where —Av = f on R^ with \v\ — > as |x| — > 00. 
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3 The small parameter 77, degeneration of the H ^-norm, and 



the rescaling of (1.2) 



We introduce a new parameter r] controlling the vanishing volume. That is, we consider the 
total mass to be rf^M, for some fixed M, and rescale as 

u 

This will facilitate the convergence to Dirac delta measures of total mass M and will lead 
to functionals defined over functions Vrj : T" — > {0, l/r/"}. Note that this transforms the 
characteristic function u of mass r]"'M to a function ti.^ with mass M, i.e., 



u = rfM while / = M. 

On the other hand, throughout our analysis with functions taking on two values {0, 1/r/"}, 
we will often need to rescale back to characteristic functions in a way such that the mass is 
conserved. To this end, let us fix some notation which we will use throughout the sequel. 
Consider a collection : T" — > {0, 1/??"} of components of the form 



< = ^Xa,, (3.1) 



i 

where the Ai are disjoint, connected subsets of T". Moreover, we will always be able to 
assum^ without loss of generality that the Ai have a diameteij^ less than 1/2. Thus by 
associating the torus T" with [—1/2, 1/2]", we may assume that the Ai do not intersect the 
boundary 3[— 1/2, 1/2]" and hence we may trivially extend to R" by defining it to be zero 
for X ^ Ai. In this extension the total variation of calculated on the torus is preserved 
when calculated over all of R". We may then transform the components v*, to functions 
: R" — > R by a mass-conservative rescaling that maps their amplitude to 1, i.e., set 

z;(x):=r/X(r?x). (3.2) 



We first consider the case n = 3. Consider a sequence of functions of the form (3.1). 
The norm \\v^ — J can be split up as 



oo 



'Jr, - i V'nfH-i{T3) = ^^vli{x)vl^{y)Grj^3{x - y)dxdy 

i=l 



CxD „ „ 

^ y^3y^3^^(^)^^(y)'^T3(a:-y)da;dy. (3.3) 



*We will show in the course of the proofs that this basic Ansatz of separated connected sets of small 



diameter is in fact generic for a sequence of bounded mass and energy (c/. Lemma 5.2 I. 
^For the definition of diameter, we first note that the torus T" has an induced metric 

d{x, y) := min{|a; — y — k\ : fc G Z"} for x,y £ T". 

The diameter of a set is then defined in the usual way, 

diam A :— sup{d(a;, y) : x,y £ A\. 
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As we shall see (c/. the proof of Theorem 4.2 ), in the limit ^ it is the first sum, containing 



the diagonal terms, that dominates. For these terms we have 
Ibn - f <IIh-i(t3) = / / vU^)vUy) G {x - y)dxdy 

vl^{x)vl^{y) -^\x - y\~^ dxdy + I I vl^{x)vl^{y) g^^\x - y) dxdy 

7] " I I zUx/r])zUy/ri)^\x-y\~^dxdy + 
I-rsJr3 47r 
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+ / / v^{xy {y) g^^\x - y) dxdy 



vli{x)vli{y) g^'^\x - y) dxdy 



Air 

(3) 



V ^\\zv\\h~^r3) + J^^ J^^^vli{x)vli{y) g^^>{x - y) dxdy. 



(3.4) 



This calculation shows that if the transformed components converge in a 'reasonable' sense, 
then the dominant behavior of the H~^-novm of the original sequence v is given by the term 



(R3) 



o 



This argument shows how in the leading-order term only information about the local behavior 
of each of the separate components enters. The position information is lost, at this level; we 
will recover this in the study of the next level of approximation. 
Turning to the energy, we calculate 



/f-i(T3) 



■ '"I|//-1(T3) 
|2 



£{u) = / |Vu| + 7 — / 

= rf" I \Vv\ + ■-frf \\v - j 

Note that if Vr, consists of = 0(1) particles of typical size 0(r/), then 

T] / \Vv\ ~ 0(1). 
Jt3 



(3.5) 



Prompted by (3.4), we expect to make both terms in (3.5) of the same order by setting 



Therefore we define 



1 



E^^(0 := -^£{u) 



1 



T] /t3 |V^;| +v\\v-f ^^IIh-i(t3) if ^ e W(T3; {0, l/r/^}) 
oo otherwise. 
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We now switch to the case n = 2. Here the critical scahng of the H ^ in two dimensions 
causes a different behavior: 



vUx)v'Jy) Gt2(x - y) dxdy 
1 



, , vU x)vU y) \og\x -y\ dxdy + I / vUx)vUy) g^'^\x - y) dxdy 
27r ./t2 Jt2 Jt2 Jt2 



1 

1 f~f A\ 1 



4(a:)4(y)log|??(a; - y)| + / / vUx)vUy) g^'^\x - y) dxdy 
R2 JR2 Jt2 Jt2 



o .1 4 ^ogrj-— / log dxdy 

V -/r2 7 ^R2 



+ J^^ J^^vij{xy^{y) g^^\x - y) dxdy 

= ^(/ 4) l^^SVl-T^-f I Zn{x) z'Jy) log \x -y\ dxdy 
V^R2 7 ^R2 Jr2 

+ vijix)vl^iy) 9^^\x - y) dxdy. (3.6) 

By this calculation we expect that the dominant behavior of the H'^-noivn of the original 
sequence v is given by the term 

Note how, in contrast to the three-dimensional case, only the distribution of the mass of v 
over the different components enters in the limit behavior. Note also that the critical scaling 
here is |log ri\. 

Following the same line as for the three-dimensional case, and setting 



we calculate 



v = ^, (3. 



"||//-1(T2) 



£{u) = / |Vn| -|- ^\\u — j-\ 

j^^ \Vv\ + 7^^II^^-/^^IIh-i(t2) 



T] r] 



Following (3.6), (3.7), in order to capture a nontrivial limit we must choose 

1 



^ |logr/| r/^ 



With this choice of 7, we define 



P2d, , ^ /^/T2|V^^| + |logr?rMl^^-f^lli-i(T2) ift^G W(T2;{0,1/??2}) 

E {v) := -S{u) = < ■ 

V 00 otherwise. 
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4 Statement of the main results in three dimensions 

We now state precisely the r-convergence results for E^"^ in three dimensions. Both our F- 
limits will be defined over countable sums of weighted Dirac delta measures Yli^i fn^^x^- We 
start with the first-order limit. To this end, let us introduce the function 

ef{m) ■= inf |Vz| + |k||i-i(R3) : z G BV{B?; {0, 1}), z = m| . (4.1) 

We also define the limit functional 



X^^x Cq'^C™'*) a V = X^^]^ fn^^x^^ {^*} distinct, and m* > 
oo otherwise. 



Remark 4.1. Under weak convergence multiple point masses may join to form a single 
point mass. The functional Eq^ is lower- semicontinuous under such a change if and only if 
the function Cq^ satisfies the related inequality 



oo oo 



mM < Y.^l^'i^')- (4-2) 



i=l 1=1 

The function Cq^ does satisfy this property, as can be recognized by taking approximating 
functions and translating them far from each other; the sum - is admissible and its 
limiting energy, in the limit of large separation, is the sum of the individual energies. 

Having introduced the limit functional Eq"^, we are now in a position to state the first 
main result of this paper. 

Theorem 4.2. Within the space X, we have 

£f as ry ^ 0. 

That is, 



(Condition 1 - the lower hound and compactness) Let Vrj be a sequence such that the 

i3d 

is countable, and 

I 



sequence of energies E^^(t;^) is bounded. Then (up to a subsequence) vq, suppvo 

\i^miEl\vr,)>El\vo). (4.3) 



^The definition of requires the point mass positions to be distinct, and tlie reader miglit wonder wfiy 
tfiis is necessary. Consider tfie following functional, which might be seen as an alternative. 



Y:Zi ^o^im') if « = E^i m'S,. with > 0, 
oo otherwise. 



This functional is actually not well defined: the function v will have many representations (of the type 5 = 
aS + {1 — a)5, for any a G (0, 1)) that will not give rise to the same value of the functional. Therefore the 

functional Eq'' is a functional of the representation, not of the limit measure v. The restriction to distinct a;' 
eliminates this dependence on representation. 
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(Condition 2 - the upper hound) Let E^{vq) < oo. Then there exists a sequence Vrj 
such that 



limsupE^''(t;^) < Eo''(t;o). 

Note that the compactness condition which usuahy accompanies a Gamma-convergence 
result has been built into Condition 1 (the lower bound). The fact that sequences with 
bounded energy E^^ converge to a collection of delta functions is partly so by construction: 
the functions are positive, have uniformly bounded mass, and only take values either or 

Since — > 0, the size of the support of shrinks to zero, and along a subsequence 
converges in the sense of measures to a limit measure; in line with the discussion above, this 



i3d/ 



limit measure is shown to be a sum of Dirac delta measures (Lemma 5.1). 



We have the following properties of Cg'^, and a characterization of minimizers of Eq^. The 



proof is presented in Section 5.4 



Lemma 4.3. 1. For every a > 0, Cq^' is non-negative and bounded from above on [a, cxd). 

2. is strictly concave on [0,27r]. 

3. If {m*}jgN with ^ • m* < oo satisfies 



m 



(4.4) 



i=l 



then only a finite number of m* are non-zero. 



Note that the limit functional Eq"^ is blind to positional information: the value of Eq'^ 
is independent of the positions of the point masses. In order to capture this positional 
information, we consider the next level of approximation, by subtracting the minimum of Eg'^ 
and renormalizing the result. To this end, note that among all measures of mass M, the 
global minimizer of Eq^ is given by 



mm 



i3d/ 



T3 



M} =ef{M). 



We recover the next term in the expansion as the limit of E: 
that is of the functional 



3d 

V 



Cq^, appropriately rescaled. 



^0 



T3 



If this second-order energy remains bounded in the limit rj 
vq = ni^S^i necessarily has two properties: 



0, then the limiting object 



1. The limiting mass weights {m*} satisfy (4.4); 



2. For each m*, the minimization problem defining eg'^(m*) has a minimizer. 

The first property above arises from the condition that E^^(v^) converges to its minimal value 
as r/ — > 0. The second is slightly more subtle, and can be understood by the following formal 
scaling argument. 
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In the course of the proof we construct truncated versions of Vrj, called , each of which is 



localized around the corresponding limiting point and rescaled as in (3.2) to a function z*. 



For each i the sequence is a minimizing sequence for the minimization problem Bq {m^), 
and the scaling of F^*^ implies that the energy E^^(t;^) converges to the limiting value at a rate 
of at least 0{r]). In addition, since converges to a delta function, the typical spatial extent 
of supp is of order o(l), and therefore the spatial extent of supp is of order 0(1/7?) . If the 
sequence does not converge, however, then it splits up into separate parts; the interaction 
between these parts is penalized by the H~^-noim. at the rate of 1/d, where d is the distance 
between the separating parts. Since d = o(l/ry), the energy penalty associated with separation 
scales larger than 0{r]), which contradicts the convergence rate mentioned above. 

This is no coincidence; the scaling of F^'^ has been chosen just so that the interaction 
between objects that are separated by 0(l)-distances in the original variable x contributes an 
0(1) amount to this second-level energy. If they are asymptotically closer, then the interaction 
blows up. 

Motivated by these remarks we define the set of admissible limit sequences 
Ai := llm-^jjgN : > 0, satisfying (4.4), such that eQ'^(m*) admits a minimizer for each z| 
The limiting energy functional Fq^ can already be recognized in the decomposition given 



p3d 



by (|3.3|) and (|3.4|). We show in the proof in Section |5| that the interfacial term in the energy 

as is the highest-order term in 



is completely cancelled by the corresponding term in Cq 



3d 



the expansion of — ^ ' 
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What remains is a combination of cross terms. 



t3 Jt 



vlj{x)v^^{y)GT^3{x - y) dxdy, 



and lower-order self-interaction parts of the H ^-norm. 



1=1 



I^J^^^vl,ix)vl^{y)g^^'ix - y) dxdy. 



With these remarks we define 

- 00 

n 

Gx3(x* — x^) ii V = ^ m^5^i with {x*} distinct, {m*} G M 



i=l 



00 



1=1 

otherwise. 



We have: 



Theorem 4.4. Within the space X, we have 



-3d r 



-3d 



as rj 



0. 



That is, Conditions 1 and 2 of Theorem 4.2 hold with and Eq replaced with F^ and Fq 
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The interesting aspects of this hmit functional Fq*^ are 

• In contrast to Eq^, the functional Fq'^ is only finite on finite collections of point masses, 



which in addition satisfy two constraints: the collection should satisfy (4.4), and each 



weight should be such that the corresponding minimization problem (4.1 ) is achieved. 
In Section [7] we discuss these properties further. 

The main component of Fg^ is the two-point interaction energy 

id- i+i 

This two-point interaction energy is known as a Coulomb interaction energy, by reference 
to electrostatics. A similar limit functional also appeared in 



5 Proofs of Theorems 14.21 and 14.4 



5.1 Concentration into point measures 

Lemma 5.1 (Compactness). Let he a sequence in BV{T^; {0, I/t]^}) such that both J^ps v-q 
and E^'^(f^) are uniformly bounded. Then there exists a subsequence such that v^^ vq as 
measures, where 

oo 

VQ■=^m'■5^^, (5.1) 

i=l 

with m* > and x* G T'^ distinct. 

Note that we often write "a sequence v^" instead of "a sequence — > and a sequence Vn^ 
whenever this does not lead to confusion. The essential tool to prove convergence to delta 
measures is the Second Concentration Compactness Lemma of Lions [18] . 



Proof. The functions Wn ■= rjVri satisfy — > in L^{T'^), and |Vw^| = r] \Vvri\ bounded in 
L^(T^). On the other hand, by definition, one has wfj^"^ = which is bounded in L^(T^). 
Hence we extract a subsequence such that Vn vq as measures. Lemma I.l (i) of [IB] (with 



m = p = l,q = 3/2) then implies that vo has the structure (5.1 ). □ 



The proof of the two lower-bound inequalities uses a partition of supp Vr^ into disjoint sets 
with positive pairwise distance. This division implies the inequality 




and is an important step towards the separation of local and global effects in the functionals. 
The following lemma provides this partition into disjoint particles. 

Lemma 5.2. Continue under the conditions of the previous lemma. For the purpose of 
proving a lower bound on E^'^(u^) and f'^^{vrj) we can assume without loss of generality that 
for some n G N 

n 
i=l 
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with w-liminf ^_^o ^ t^q^x^ o,s measures, dist(suppt'* ,suppw^) > for all i ^ j, and 
diamsuppf* < 1/4. In addition, for the lower bound on F^'^(f^) we can assume that for each 



i, iT^o^x^ ^''^d that there exist ^* G and a constant C* > such that 



x-Ci^\'^vij{x)dx < C'r] 



T3 



(5.2) 



The proof of Lemma 5.2 is given in detail in Section 5.4 A central ingredient is the 
following truncation lemma. Here $7 is either the torus T'^ or an open bounded subset of R^. 



Lemma 5.3. Let n ^ISi be fixed, let oo, and let Uk S BV{0,; {0, a^}) satisfy 



\Vuk\ = o(afc), 



(5.3) 



and converges weakly in X to a weighted sum 



1=1 



where m' > and the G 17 are distinct. Then there exist components u\ £ BV{^1; {0, Ofc}), 
i = 1 . . .n, satisfying diam supp u^, < 1/4, inf^ infj^j dist(supp u^, supp u-j^) > 0, and 



w-liminf li^ > m^6j.i, 

fc— >oo 

in the sense of distributions. In addition, the modified sequence = satisfies 

1. Uk < Uk for all k; 

2. limsupfc^o^ Jiuk-Uk) < 

3. There exists a constant C = C{n) > such that for all k 

\Vuk\ < J \Vuk\-C\\uk-Uk\\L3/^{n)■ 



{5.A) 



(5.5) 



The essential aspects of this lemma are the construction of a new sequence which again 



lies in BV{^; {0,0^}), and the quantitative inequality (5.5) relating the perimeters. 
5.2 Proof of Theorem lOl 

Proof. (Lower bound) Let Vrj be a sequence such that the sequences of energies E^^{vri) 
and masses /r^a are bounded. By Lemma 5.1 a subsequence converges to a limit vq of the 
form (5.1 ). By Lemma 5.2 it is sufficient to consider a sequence (again called Vfj) such that = 
X^iLi ■yj; with w-liminf^^o ^^J; ^ ^o^x^^ suppuj^ C B{x^,l/A), and dist(suppu,^,suppt>^) > 
for all i ^ j. Then, writing 

we have 



4(y) ■.= V^vUx' + riy), 



(5.6) 



iT3 Jr^ 



and 



/r3 
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and by K^ 



h-^t3)=V 'II^^II//-i{r3) + j^^^ j^^vl^{x)vl^{y)g^^\x -y)dxdy. 



For future use we introduce the shorthand 



/t3 Jrs 



Then 



n n „ „ 

^=l ^,j=l -^T^ iT3 



dx dy 



E 



\Vzl\ + \\z. 



i ||2 



R3 



i=l 



j ^vl^[x)vli[y)g^^\x-y)dxdy + riY^ j ^ j y^(x)vl^{y)GT.3{x - y) dx dy 



> 



^ ef (m;) + inf 5^^^ ("i;) % inf Gt3 ^ . 



(5.7) 



i=l 



i=l 



Since the last two terms vanish in the hmit, the continuity and monotonicity of e^^ (a conse- 



quence of Lemma 4.3) imply that 



hmjnf E^^K) >Y,el' limjnf mM > ^e^^^K) > EI'{vo). 
i=i \ / j^^i 



(Upper bound) Let vq satisfy Eq^(uo) < oo. It is sufficient to prove the statement for 
finite sums 



i=l 



since an infinite sum vq = Yl"^! ^x^- can trivially be approximated by finite sums, and in 
that case 

/ n \ n oo 



i=l 



1=1 



To construct the appropriate sequence Vr^ ^ vq, let e > and let be near-optimal in 
the definition of eQ*^(m*), i.e., 



\Vz'\ + \\z 



i\\2 



R3 



lH-i{R3) 



n 



(5. 



By an argument based on the isoperimetric inequality we can assume that the support of 
is bounded. We then set 

-^Ji^-if^ -^^^ (5.9) 



vUx):=r] '^z\r] \x - x")) 
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so that 



T3 



= m . 



Since the diameters of the supports of the tend to zero, and since the are distinct, 
Vr^ := is admissible for E ^^ w hen rj is sufficiently small. 



Following the argument of (5.7), we have 

n 



i=l 



i\\2 



R3 



i=l '^'^ '^^ i,7=l '^^ 



and thus 



limsupE^''(t;^) < E^''(t;o) + e. 



The result follows by letting e tend to zero. 



□ 



5.3 Proof of Theorem 14.41 

Proof. (Lower bound) Let Vr) = Yl^=i'^rj be a sequence with bounded energy F^'^(t>^) as 
given by Lemma |5.2[ converging to a vq of the form 



i=l 



where ttiq > and the are distinct. Again we use the rescaling (5.6) and we set 



R3 



Following the second line of (5.7) we have 



^3d 



Vr 



n 



i=l 



R3 



iVz* I + \\z. 



T3 
i ||2 



o3d / i - 



+ 



.i=l 



,1=1 



i=l T3 i.j=l 



(5.10) 



Since the first two terms are both non-negative, the boundedness of ^^{vn) and continuity 
of Cq"^ imply that 



^<Y,ef{m^)-efr£m'] = lim 
i=i \i=i J 



<o, 
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and therefore the sequence {m*} satisfies (4.4). 

By the condition (5.2) the sequence 2* is tight, and since it is bounded in BV{Ii^; {0, 1}), 
a subsequence converges in L^(R^) to a hmit (see for instance Corollary IV. 26 of f8]). We 
then have 



0< / |Vz^| + 114111-1(1,3) -eg^M 



< liminf 



R3 



lim Cn^ (m. 



V"-ri) 



which implies that Zq is a minimizer for eQ'^(m*). 
Finally we conclude that 



+ 5^ / / vUx)vUy)GT^z[x - y)dxdy\ 

n n 



i=l 



(Upper bound) Let 



vo 



i=l 



with the x"^ distinct and {m*} G 7W. By the definition of M we may choose 2* that achieve 
the minimum in the minimization problem defining e^{m!')] by an argument based on the 
isoperimetric inequality the support of is bounded. 

Setting vl^ by (5.9), for 77 sufficiently small the function v^q := X^iLi admissible for 
F^^, and v„ 



vq. Then following the second line of (5.7), we have 



limF^^(^;,) = F^^(^;o). 



□ 



5.4 Proofs of Lemmas 15.21 and 15.31 



For the proof of Lemma 5.2 we first state and prove two lemmas. Throughout this section, if 
i? is a ball in and A > 0, then XB is the ball in obtained by multiplying i? by A with 
respect to the center oi B] B and XB therefore have the same center. 

Lemma 5.4. Let w G BV{Br\ {0, 1}). Choose < r < R, and set A := Bji \ Bj.. Then for 
any r < p < R we have 



n'^idBpHsuppw) 1 



/ |Vtt;| + + w 

J A J A 
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Proof. Let P be the projection of onto Br- For any closed set D C with finite 
perimeter, the projected set P(A f] D) is included in Ei, U E^^ where the two sets are: 

• The projected boundary Ei, := P{A n dD); since P is a contraction, n'^{Ei,) < n^{An 
dD); 

• The set of projections of full radii Ej. := {x £ dBr : Xx G D for all 1 < A < R/r}, for 
which 

n\Er) = ^§^^'({Ax :xeEr,l<X< R/r}) < ^^^^C^D n A). 

Applying these estimates to D = supp w we find 

n^idBpHsuppw) _ n^{P{dBpnsuppw)) 
W{dBp) ~ W{dBr) 

n^{p{Ansuppw)) 



W^idBr) 

1 [av2,.^^ ^ , n-'idBr,^^, 



< 

- W{dB) |^'(^n^suppti;) + '^^3^^|' /:^(^nsuppu>: 

and this last expression implies the assertion. □ 

Lemma 5.5. There exists < a < 1 with the following property. For any w G BV{Bfi] {0, 1}) 
with 

^7277^/ _|Vu;| + / _w<\, (5.11) 
there exists a < P < 1 such that 

2\\dBpR\\{suppw) < \Vw\. (5.12) 

J Bji\Bi3ji 

Proof. By approximating (see for example Theorem 3.42 of |1 ) and scaling we can assume 
that w has smooth support and that R = 1. Set < a < 1 to be such that 

^l-^ = n'{dB^), (5.13) 

where C is the constant in the relative isoperimetric inequality on the sphere 5^: 
mm{n'^{D n 5^), 7^2(52 ^ < c{Ti}{dD n S'^)f . 



We note that the combination of the assumption (5.11) and Lemma 5.4 implies that when 
applying this inequality to D = suppw, with S"^ replaced by dBis, the minimum is attained 
by the first argument, i.e. we have 

n'^{D n dBi^s) < c{n\dD n dBi^s)f- 

We now assume that the assertion of the Lemma is false, i.e. that for all a < r < 1 

< 2\\dBr\\{D) - \\dD\\{Bi\Br). (5.14) 
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Setting f{s) := H^dD n dBi^s) we have 



f f{a)da= [ n\dD n dBr) dr < [ _ \Vw\^ 2\\dBi^s\\{D). (5.15) 

Jo Jl-s Jbi\Bi^s 



f{a) da 

By the relative isoperimetric inequality we find 



/ /(a) da < 2\\dBi^s\\{D) < 2C{n\dD n dBis))^ = 2Cf{sf. 
Jo 

Note that this inequality implies that / is strictly positive for all s. Solving this inequality 
for positive functions / we find 



l-a 



f{a) da > 



SC 



2n\dB^) > 2\\dBj{D) ^ 



l-a 



f{a) da, 



a contradiction. Therefore there exists an r =: (3R satisfying (5.12), and the result follows as 
remarked above. □ 



Proof of Lemma 5.3: Let a be as in Lemma 5.5 Choose n balls i?*, of radius less than 
1/8, centered at {x*}"^;^, and such that the family {25*} is disjoint. Set Wk ■= a^^Uk-, and 
note that for each z, 

C 



_Wk < 



/ IVUfcl + / Uk \ , 

Jn Jn ) 



and this number tends to zero by (5.3), implying that the function Wk on B^ is admissible for 
Lemma 



5.5 



5.5 



For each i and each A;, let be given by Lemma 

2||a/?ii?*||(suppnfc) <a,i||Vt.fc||(i?*\45^). 



so that 



(5.16) 



Now set u\ := UkXpi^B^- ^^"^ '■= Yll=i ^k- Then for any open A G ^ such that G A, 
liminf I u\ = liminf Uk > liminf Uk > / m^6r^j{An aB^) > m*, 

k—>oo J k—>oo J k—>oo J 7— ' 

which proves (5.4); property [2] follows from this by remarking that 

limsup / {uk — Uk) = li™ Uk — liminf Uk m-' — > 

fc^oo Jn k~^°°Jn fc-*oo Jn ~[ ~[ 

The uniform separation of the supports is guaranteed by the condition that the family {2B'} 
is disjoint, and property [T] follows by construction; it only remains to prove (5.5). 
For this we calculate 

„ n n 

/ iVufcl = \\Vuk\\([j PIB'^ + akJ2\\9PiB'\\{suppuk) 
•^^ i=i i=i 

■^^ i=i i=i 

\Vnk\--\\Vuk\\{n\[jf3iB') 

i=l 

< I iVufcl - Cfcllufc -/^^Ufc||^3/2(^^) (5.17) 
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Here the constant Ck is the constant in the Sobolev inequality on the domain Aj. := r2\Uj/3^i?*, 

The number Ck > depends on k through the geometry of the domain Ak- Note that the 
size of the holes is bounded from above by and from below by aB''. Consequently, 

for each ki and /c2 there exists a smooth diffeomorphism mapping A^-^ into ^^j, and the first 
and second derivatives of this mapping are bounded uniformly in ki and k2- Therefore we 
can replace in (5.17) the /c-dependent constant by a ^-independent (but n-dependent) 
constant C > 0. 

Note that since Uk is bounded in L^, 



-3/2, 



Uk 



|3/2 _ -111 II 







as k ^ oo. 



(5.18) 



Continuing from (5.17) we then estimate by the inverse triangle inequality 
\Vuk\ < iVufcl - C||nfc||^3/2(^^) + 1^^1/3 ll^fcllLi(Afe) 



/ iVufcl - Cllufcl 
Jo. 

/ iVufcl - C\\uk\ 
Jo, 



L3/2{Afe) 



+ 



c 



|3/2 



L3/2(A,) 



/3 l/2ll-'=ll.3/2(^,) 



1 1/2 



s y |Vufc| - c"||tifc||^3/2(^^) = y \vuk\-c'\ 



Uk - Uk\\L3/2^Qy 



This proves the inequality (5.5). 



□ 



Proof of Lemma 5.2 By Lemma 5.1 and by passing to a subsequence we can assume 
that v-q converges as measures to vq. We first concentrate on the lower bound for E^^. 

Fix n G N for the moment. We apply Lemma 5.3 to the sequence Vrj and find a collection 
of components f J^, z = 1 . . . n, and Vjj = u^, such that 



w-limmf v„ > y m 6^i, 
i=i 



and 



T3 



C\\n 



T3 



ri ~ f^»?llL3/2(T3)- 



Setting := Vrj — Vrj we also have 



I ~ X ~ l|2 
ll'j? - f Vrj\\H-i 



{T3) 



Vr^{x)vrj{y)GT.-i{x - y) dxdy 



T3 JT3 



?'»?(a^)^^r;(y)G'T3(x - y) dxdy 



rri{x)r.ri{y)GT^3{x - y) dxdy 
— ~ / ^?yllii'-i(T3) ~ 2inf Gx3||r,,||ii(T3)||w^||Li(T3 
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Therefore 



Assuming the lower bound has been proved for Vri, we then find 



(5.19) 



liminf (ii^) > liminf (ti^) — C'||r^||/^i(x3) 



> E^^ ( w-liminf Vr,) - C hm / v„ + C" hm inf 



1=1 



m . 



i=n+l 



Taking the supremum over n the lower bound inequality for follows. 
Turning to a lower bound for F^*^, we remark that by Lemma 



with non-zero weight m* is finite. Choosing n equal to this num 



4.3 



the number of in (5.1 ) 



aer, we have 



lim 



> liminf 



n „ n 

> hm inf / v'L>'S^ 

~ ^ '7-0 jr^Z ~ ^ 

i = l 1 = 1 



m 



lim 



and therefore v]^ m^Sr^i^ and Jr^s 0. Then 



( |5.19| 1 

ri 



C|k»?llL3/2(T3) + — lkr?llLi(T3) 



+ 



c 



^0 



T3 

2/3 ^ L + C 



^0 



T3 



T3 



(5.20) 



Here L is an upper bound for eQ*^' on the set [inf J VrpOo) (see Lemma 4.3) and in the 
passage to the last inequality we used the triangle inequality for || • ||^3/2 and the fact that by 
construction, r^j takes on only two values. For sufficiently small r], the last two terms add up 
to a negative value, and therefore we again have F^^{vri) < F^^{vri). Because of the choice of 
n we have Vrj vq; if we assume, in the same way as above, that the lower bound has been 
proved for Vr], we then find that 



liminf F^^(t;^) > liminf F^^(i)^) > Ff (z;o). 
t;— >0 ' r|—^0 ' 
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For use below we note that 
1 



V 



i=l 



i= 

n 



1 



+ ||wJ,||^-i(t3) - -co 



|Vi;;| + III;; 



3d 



X3 JT3 



vl,{x)vUy)GT^s{x - y) dxdy 



(5.21) 



In the calculation above, and in the remainder of the proof, we switch to considering f * defined 
on instead of T'^. Since the terms in the first sum above are non-negative, boundedness 
of f^^{vr]) as t] ^ implies the boundedness of each of the terms in the sum independently. 
We now show that when F^'^(t>^) is bounded, then for each i 



G : 



^^|^w^(x)(ix = 0(7/^) as 7? ^ 0. 



R3 



(5.22) 



Suppose that this is not the case for some i; fix this i. We choose for the bary center of u^, 
i.e. 



/ xvl. 


(x) dx 






f 









(5.23) 



Since we assume the negation of (|5.22|), we find that 

\x — iri\^ vl^{x) dx » ?7^. 



R3 



Note that by (5.23) and the fact that ^ x*, 

lim Pn = 0. 

Now rescale by defining ^^^{x) := p^vl^{£,r) + p-qx). The sequence C-n satisfies 
1. C,G W(R3,{0,/53r/-3}); 



(5.24) 



(5.25) 



2. 



3. 



4. 



R3 



R3 



Pv V^R3 



VCt?! + llCr,|lH-i(R3) I^^^l +^II^^^IIh-i(r3)> and 



R3 



/R3 



l^l^ Crj{x) dx = 1. 
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The first three properties imply that the sequence is of the same type as the sequence Vrj 
in the rest of this paper, provided one replaces the small parameter rj by the small parameter 
fj := rj /p-n- The fourth property implies that the sequence is tight. By the third property 
above, (5.21), and (5.25), the boundedness of F^^ translates into the vanishing of the analogous 
expression for (ri- 



limsupi/ |VC,| + ||C.II?,-i(R3)-^er(/ 



0. 



(5.26) 



We now construct a contradiction with this limiting behavior, and therefore prove (5.22). 



Following the same arguments as for Vrj we apply the concentration-compactness lemma 
of Lions [T8] to find that the sequence (^^ converges to (yet another) weighted sum of delta 
functions 

oo 

/i := ''^m^Syj with Ai(R^) = rn\ 

where m-' > and G are distinct. Since 

/xd.W^lim/' .C,W<ix = and Z H= d.W = lim / |.f C, W <i. = 1. 

at least two different m-' are non-zero; we assume those to be j = 1 and j = 2. 

We will need to show that the number of non-zero m-^ is finite. Assuming the opposite for 
the moment, choose n G N so large that 

n 

j;ef(m^)>e3V'); 
i=i 

this is possible since the re exist no minimi zers for eQ^{m^) with infinitely many non-zero 
components (Lemma 4.3). We apply Lemma 5.3 to find a new sequence Cr; = Z]j=i Cn^ where 
C^^m^6yj. Then 

n „ 

>liminf^eg'i(/ C^) > (m^), 



which contradicts (5.26); therefore the number of non-zero components m-' is finite, and we 
can choose n such that m* = YTj=\ I iCv ~ Cr;) — *■ 0. 
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To conclude the proof we now note that 



^|V<.,|+||(,||l,_,„,-'^e«(/^_C,; 

n 

- Xl{ / l^^^l + II^^IIh-1{R3)} +2(C,),C^)h-1{R3) + C\\Cr^ - C,,||l3/2(r3) 
7 = 1 "'R-^ 



> 



R3 ^ VJr3 



+ 



+ 2(C?,5 Cr,)H-i(R3) + C'lIC?? - C??IIl3/2(R3) 



> 



~ \2/3 1 

r/ \Jr,3 ' ' / 27r yR,3 7^3 \x - y\ 



dxdy. 



Since hm^^o J-RziCri — Cri) = 0, the first two terms in the last line above eventually become 
positive; the final term converges to (27r)~^m"'^m^|y^ — y'^\^^ > 0. This contradicts (5.26). □ 



5.5 Proof of Lemma 14.31 

Let Zn be a minimizing sequence for eQ^{m). The functions 

\(1 + e/m)^''^ 

are admissible for eQ'^{m + e) for all e > —m. Since the functions 
Us) := / |V4| + \\z: 

Jrs 

satisfy 



2nll//-i(R3) 



(1 + e/m)2/3 / |Vz„| + (1 + e/mf/^z„ 



fn{e) = /„(0) + + ^lkn||H-i(R3) 



uniformly in n, we have for all e > 0, 

el''{m + e) < inf < (m) + ^ef (m)^ + ^e-j^^M f^)' + of 



We deduce that 



e3d(^ + e) - e^'iCm) < ^el'^{m) e + 0(e2). 



(5.28) 
(5.29) 



By (|4.2|), we find that for any m > 1 and any positive integer n, we have 

' m — l' 



3d/ 



o3d 



n 
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By taking n such that ^ G [1,2], we have 

el'^im) < eg^(l) + Cn 

where C denotes a uniform bound for Cq^ on the interval [1,2]. By cho ice of n we have for 
some constant C", eg'^(m) < eo^(l) + Cm. Combining this with (5.29), we find that Cq*^ is 
bounded from above on sets of the form [a, oo) with a > 0. 

For the concaveness of e^'^, note that under a constant-mass constraint / |Vz| is minimal 
for balls and ||^||_h--i(r3) is maximal for balls (see e.g. ^Oj for the latter). Setting m = J z and 
= 3m/47r, we therefore have 



f 10 2 2 /■ 10 2 



and an explicit calculation shows that the right-hand side is negative iff m < 27r. From (5.27) 
we therefore have for all m < 27r and all e > —m, 

el'^{m + e) < el'^{m) + inf a„e - be'^ + ce^ 

n L 

where a„ is a sequence of real numbers, and 6, c > 0. Writing this as 

eo'^(m) < inf \ 60*^(771,0) -|- inf a„(r77 — 7779) — 6(777 — mo)^ + c(m — 7770)^ \ , 

moe(0,27r) I n I J J 

we note that for each 7779 the expression in braces is strictly concave in m for 1 777 — 7779 1 < b/3c; 
since the infimum of a set of concave functions is concave, it follows that the right-hand side 
is a concave function of m. Since equality holds for 7779 = 777, eQ*^ is therefore concave for 
777 < 27r, and eo*^"(77i) < for 777 < 2tt. 

Finally, part 3 follows from remarking that if (say) m^,m? £ (0,27r), then 

^ [elHm' +e) + erim' - e)) [^^ = ef'{m') + ef'{m') < 0. 

Therefore the sequence (777^, 777^^, . . .) is not optimal, a contradiction. It follows that there 
can be at most one 777* in the region (0,27r), and since the total sum is finite, the number of 
non-zero 777* is finite. □ 



6 Two dimensions 

All differences between the two- and three-dimensional case arise from a single fact: the 
scaling of the is critical in two dimensions, making the two-dimensional case special. 

6.1 Leading-order convergence 

The first difference is encountered in the leading-order limiting behavior. As we discussed in 
Section [sj the leading-order contribution to the i^~^-norm involves the masses of the particles 



instead of their localized H ^-norm (see (3.7)). For the local problem in two dimensions we 
therefore introduce the function 

t2 



eoim) := - 

2 



- + inf|/" \Vz\: z e BV(R^;{0,1}), [ z = m] (6.1) 



m 

h 2y'irm. 
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Note that the minimization problem in (6.1 ) is simply to minimize perimeter for a given area, 
and a disc of the appropriate area is the only solution. Thus the value of eQ^(m) can be 
determined explicitly. 



The function Cq does not satisfy the lower-semicontinuity condition (4.2) (c/. Remark 



4.1). We therefore introduce the lower-semicontinuous envelope function 



e2d(m) := inf <^ ^eg^(mJ) : > 0, ^ 



m-^ = m 



(6.2) 



The limit functional is defined in terms of this envelope function: 



Eoo 
i=l 



m 



) if i; = Yli^i ^^^x^ with {x'} distinct and m* > 



CXD 



otherwise. 

Theorem 6.1. Within the space X, we have 

as r/ — > 0. 

That is, conditions 1 and 2 of Theorem 4.2 hold with E^'^ and Eq"^ replaced by E^^ and Eq^. 



p2d r r.2d 



The proof follows along exactly the same lines as the proof of Theorem 4.2 It is in fact 
simpler, since a standard result on the approximation for sets of finite perimeter (see for 
example Theorem 3.42 of |3]) implies that, without loss of generality, we may assume that a 
sequence Vfj with bounded energy (for r/ sufficiently small) satisfies 



'Jrj = ^ f * with V, 



i=l 



1 



(6.3) 



where the sets Al^ are connected, disjoint, smooth, and with diameters which tend to zero as 
T] ^ 0. Then the following estimate holds true in two dimensions: 

oo oo „ 

j;diam(supp^;;) < r/^ ^ / \Vvi^\ < V^^'ivr,) = 0{ij), (6.4) 
1=1 i=i •'^^ 



which can be used to bypass Lemma 5.2 



6.2 Next-order behavior 

Turning to the next-order behavior, note that among all measures of mass M, the global 
minimizer of Eq*^ is given by 

m.ini^f{v):j v = =^{M). 

We recover the next term in the expansion as the limit of E^'^ — e^, appropriately rescaled, 
that is of the functional 



^f{v) := |logr/| 
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Here the situation is similar to the three-dimensional case in that for boundedness of the 
sequence F^^ the limiting weights m* should satisfy two requirements: a minimality condition 
and a compactness condition. The compactness condition is most simply written as the 
condition that 

'^{m') = el'^im') (6.5) 

and corresponds to the condition in three dimensions that there exist a minimizer of the 
minimization problem (4.1). 

In two dimensions, the minimality condition (6.6) provides a characterization that is 
stronger than the in three dimensions: 

Lemma 6.2. Let {m*}igN be a solution of the minimization problem 

min i el'^im') : > 0, ^ = M. I . (6.6) 



.i=l 1=1 



Then only a finite number of the terms are non-zero and all the non-zero terms are equal. 
In addition, if one is less than 2^^/^7r, then it is the only non-zero term. 



The proof is presented in Section 6.3 We will also need the following corollary on the 
stability of Eq'^ under perturbation of mass: 

Corollary 6.3. The function Cq*^ is Lipschitz continuous on [6, 1/6] for any < 6 < 1. 



The limit as r/ ^ of the functional F^*^ has one additional term in comparison to the 
three-dimensional case, which arises from the second term in (3.6), 



^ oo 

-y 



\x - y\ dxdy. 



(6.7) 



To motivate the limit of this term, recall that appears in the minimization problem (6.1 ), 
which has only balls as solutions. Assumin g -z^ to be a characteristic function of a ball of 
mass m', we calculate that the first term in (6.7) has the value fo{m}), where 



„ , , m / , m 
fo{m) := — 3-21og- 
ovr \ vr 

We therefore define the intended T-limit Fg"^ of F^^ as follows. First let us introduce some 
notation: for n G N and m > the sequence n (8> m is defined by 

I m 1 < i < n 

[n (g) m) := < 

n -I- 1 < i < oo. 



Let M be the set of optimal sequences for the problem (6.6): 



Ai := |n ® m : n®ra minimizes (6.6) for M = nm, and e^{m) = eo'^(m)| 
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Then define 



'n{/o(m) + m2 5(2)(o)| 



+ 



m 
~2 



V = m , {x*} distinct, n®m£M., 



OO 



(6.8) 



otlierwise. 



Theorem 6.4. Within the space X, we have 



-2d r 



-2d 



as rj 



0. 



That is, Conditions 1 and 2 o/ Theorem 6.1 /loW wit/i E^*^ anc? Eq^ replaced with ¥^ and Fq"^ 
respectively. 



The proof of this theorem again closely follows that of Theorem 4.4, The compactness 



property (6.5) in the lower bound follows by a simpler argument than in three dimensions, 



however. Using the division into components with connected support (6.3), we have 



|logr?| 
|logr?| ^ 



R2 

oo 



R2 



R2 



+ |logr/|j; e2d(/ 4 

oo 

E 



,2d 



R2 



+ |log??| 



"(/R2^y-^°'a2^^ 



(6.9) 
(6.10) 



+ X] 1 ^ ^(y) k - y| c?xdy + j^^ j^^ vl^{x) vlj{y) g^'^\x - y) dx 



i=l 
oo 

E 



J 1'2 



(x) u^(y) Gt2(x - y) dx dy. 



(6.11) 



The last two lines in the development above are uniformly bounded from below. Since F^'^(t>^) 
is bounded from above, it follows that the terms in square brackets, which are non-negative, 
tend to zero. In combination with the continuity of e^^ and this implies the compactness 



property (6.5). We also remark that because the contents of the square brackets in (6.9) and 



(6.10) are zero in the limit, we find with the aid of Lemma 6.2 that the number of concentration 



points X* in the weak limit of Vr^ is finite with equal coefficient weights. Moreover, we may 
assume that there are a finite number of different components of w^, and each must converge 
to a different x*; otherwise, the last term in (6.11) would tend to oo as ry tends to 0. 
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6.3 Proofs of Lemma 6.2 and Corollary |6.3 



The proof of Lemma 6.2 contains two elements. The first e^ 



the property that Cq'^ is concave on 



0,^ 



and convex on 



ement is general, and only uses 
^ , oo^ . This property reduces 



the possibilities to a combination of (a) a finite number of equal m* in the convex region with 
possibly (b) one m* in the concave region (see jl7| Section 5.4] for a similar reasoning). The 
second part, in which possibility (b) above is excluded, depends heavily on the exact form 
of Cq*^, and is an uninspiring exercise in estimation. 



Proof of Lemma 



6.2 



For this proof only, let us abuse notation and use x,x^,y,z to 



denote positive real numbers. We note that 

el\m) = 2^'\ f [-^) with fix) = x^ + ^. 

We therefore continue with / instead of e^. Since / is concave on (O, |] and convex on 
[|, oo), the following hold true: 

• There is at most one G (O, |); for if x^x-' G (O, |), then 



(/(x* + e) + /(x^-e)) 



e=0 



/"(x*) + r(x^)<o, 



contradicting minimality. Therefore only one non-zero element is less than |, which 
also implies that the number of non-zero elements is finite. 

• The set of elements {x* : x* > is a singleton, since the function is convex on oo). 

Therefore the lemma is proved if we can show the following. Take any sequence of the 
form 

X i = 1 

y i = 2,...,n-hl (6.12) 

i>n + 2, 



with X < 1/4 < y; then this sequence can not be a solution of the minimization problem (6.6 1 . 
To this end, we first note that 



(n + 1)/ 



n 



n + l 



V 



nf{y) 



n 



n + l 



Vy + (n + 1) 



n + l 



n 



1 



If this expression is negative, then by replacing the n copies of y in (6.12 ) by n -|- 1 copies of 
ny/{n + 1) we decrease the value in (6.6). Therefore we can assume that 

\ 2/3 



\<y<yn.{n) := (n+ 1)2/3 U 



n + l 



1 



n 



We distinguish two cases. Case one: If y + x/n < ym{n), then we compare our sequence 



(6.121 with n copies oi z := y + x/n: 

/(x) + nf{y) - nf{y + x/n) = /(x) + nf{z - x/n) - nf{z) 



x'(l + 



1 



n 



2xz + \/x + n\ z 



X 

n 



n\/z 



--■ g{x,z). 
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We now show that g is strictly positive for all relevant values of x and z, i.e. for < < 1/4 
and 1/4 + 1/n < z < ym{n). 

Differentiating g{x, z)/x we find that 

ox X n 2x'^/^ V V " / 2x^/z — - 

V n 

This expression is negative: x < 1/4 implies that 

1 ^ 1 ^ 

n 2x^1'^ ' 

and by concavity of the square root function we have 



X 1 X 

n 2^z n 



so that the last two terms in (6.13) together are also negative. 
Since g{x, z)/x is decreasing in x, it is bounded from below by 

4ff(l/4, ^) = 1 (^1 + - 2z + 2 + 4n - V^) . 

The right-hand side of this expression is concave in z, and therefore bounded from below by 
the values at2: = (l-|-l/n)/4 and at z = ym{n). The first of these is 

"K^^^) ^^^^"^ (l-yiTi) >-^+2 + 2n(l-(l-l/2n)) = l. 
For the second, the expression 

4g(l/4, yra{n)) = ^(^ + ^- 2y„(n) + 2 + 4n (^^y™(n) - ^ - \/y^(n)^ 

is positive for n = 1, 2, as can be checked explicitly; for n > 3, we estimate 2^^/^ < ym{n) < 
((n + l)/2n)^/^ and therefore 



1 + ^) - 2ym{n) + 2 + 4n ym{n) - ^ - ^/yj^ 



1 I. 1 
4 



1 /n + l\^/^ 1 



> 2 — ■ — + 2 

in 



9 /n+ 1 
> - -2 ' 



2/3 



4 V 2n y 2^2-2/3 - ^ 

The right-hand side of this expression is strictly positive for all n > 3. This concludes the 
proof of case one. 

For case two we assume that y + x/n > ym^n), set 

ny + X 
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and compare the original structure with n + 1 copies of z: 



fix) + nf[y) - (n + l)/(z) = /(x) + nf ( — z - - ) - (n + l)/(z) 

n + 1 



n 

=: z). 

Note that the admissible values for z are 

n 



[z — xf + -Jx + Vny^ (n + l)z — x — (n + 



nym(n) + x 
Vmin) < z < — < ym[n). 



n + 1 n + 1 

We first restrict ourselves to n > 2, and state an intermediary lemma: 



(6.14) 



Lemma 6.5. Let n > 2. Then for all < x < 1/4 and for all z satisfying (6.14), 

h{x, z) > min{/i(0, z), h{l/4, z)}. 

Assuming this lemma for the moment, we first remark that h{0,z) > by the bound 
z > nym{n) / {n + 1) and the definition of ym- For the other case we remark that the function 



n t-^ \fn\/ {n + l)z — x — (n + 
is increasing in n for fixed z. Keeping in mind that n > 2 we therefore have 

Ml/4, z)>^{z-lf + ^ + V2^3z - I - 3V^, 

and this function is positive for all z > 2y„i{2)/3 ~ 0.51. This concludes the proof for n > 2. 
Before we prove Lemma |6.5| we first discuss the case n = 1, for which 

h{x, z) = 2{z -xf + ^/x + V2z - X - 2a/z. 

The domain of definition of z is 

^1 



:ym(l),ym(l) 



C [0.4410,0.8821]. 



The mixed derivative hzx is negative on the domain of x and z, so that 

hz{x, z) > hz{l/4, z)=4z-l+ ^ - 

^/2z-l/4 Vz 

This expression is again positive for the admissible values of z, and we find 

h{x, z) >h [x, 5ym(l)) = 2 (5ym(l) -x)"^ + y/x+ \/ym(l) " X - 2^J\ym{l). 

Similarly this expression is non-negative for all < x < 1/4, which concludes the proof for 
the case n = 1. □ 



Proof of Corollary |6.3| Fix < 5 < 2^'^/^t: and M £ [5,1/6]; by Lemma [(^2] there exist 
n,m with M = nm such that €^{^1) = ne^{m). Note that if n = 1 then m = M > S, and 
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if n > 1 then by Lemma 6.2 m > 2 ^/^vr > 6; therefore we have m > 6 and n < M/6. Since 

-0 



Cq*^ is the pointwise minimum of a cohection of functions Bq^, local Lipschitz continuity of Cq*^ 



now follows from the same property for the functions Cq*^ on the domain [5, 1 . □ 

We still owe the reader the proof of Lemma |6.5[ 
Proof of Lemma 6.5 We first show that if 4/25 < < 1/4, then hx{x, z) < 0. We estimate 
the derivative by using the bounds on z and x: 



n + 1 ^Jn 1 

hx[x, z) = -2 [z - X) = + — = 

n ^ ' 2^y{n + l)z- X 2^ 



n + 1 ^/n 5 

2n 2v^ny„(n) 4 

Note that ym is monotonically decreasing in n, and that we can estimate ym from below by 

yminf/' = (n + 1) f yiTi - l) > (n + 1) ^= i = i/^±i. 

\^ J 2 / 2 V n 



Using n > 2 we find 



< _2. 2-2/3 f!i±lV +!i±l ^ "^" + ^ 



1/3 

^ ^ " V ^ y 2n 27^ 4 V n 

The function £ is increasing on [1, cxd), and we have 

On the remaining region < x < 4/25 the second derivative h^x is negative: 

, . N + 1 \/n 1 

i^xxyXf Z) 



n 



A{{n + l)z-xf'^ 4x3/2 



1 125 

For any fixed z, therefore, the function x ^ h{x, z) takes its minimum on the boundary, 
that is in one of the two points x = and x = 1/4. Since the first derivative is non-zero 
on [4/25, 1/4], and since the second derivative is non-zero on (0,4/25], the minimum is only 
attained on the boundary. This concludes the proof of Lemma |6.2[ □ 



7 Discussion 

The results of this work provide a rigorous connection between the detailed, micro-scale model 



defined by if^ in ( 1.2 ) and the macroscopic, upscaled models given by the limiting energies Eq^, 
Fq'', Eq^, and Fg^. We now discuss some related aspects. 

Differences between tfie two- and tfiree-dimensional cases: scaling. The consequences of 
the difference between two and three dimensions in the scaling of the //^^-norm are best 
appreciated in the Green's functions in the whole space: if we replace x by r]x, then 

logr/x = logi] + logx in two dimensions, and r = — • - — r in three. 

|??X| T] \x\ 
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The difference between the additive effects in two dimensions and the multiplicative effect in 
three dimensions is responsible for the difference between the two limiting problems: 



J \\7z\ + in two dimensions, and J |Vz| + ||z||^_i in three. 

Differences between the two- and tliree-dimensional cases: local problems. Because of this 
difference in scaling, the local energy contributions e^'^ (and Cq^) and eQ*^ are necessarily dif- 
ferent, and since the two-dimensional local problem is the isoperimetric problem, its solution 
can be calculated explicitly in terms of m. For the three-dimensional local problem we can 
only conjecture on the structure of minimizers (see below). 

In addition to this, there is a difference in the handling of the lower semicontinuity in two 
and three dimensions. This comes from the fact that the definition of Bq*^ presupposes that 
the mass of z remains localized (does not escape to infinity) while the definition of e^^ does 
not. As a result, the function e^^ already has the right lower semi-continuity properties, while 
for Cq'^ we need to explicitly construct the lower-semicontinuous envelope function Cq"^. 

Some of the other differences are only apparent. For instance, the requirement, in the 
definition of Fg^, that for each m* the minimization problem eQ^{m'^) admits a minimizer, is 
mirrored in two dimensions by the compactness property eQ'^(m*) = eQ'^(m*). The reduction 



to 'blobs' of bounded and separated support (Lemma 5.2) is immediate in two dimensions, 
since it follows from the vanishing of the perimeter. 

Minimizers of the local problem in three dimensions. For the three-dimensional minimiza- 
tion problem (4.1 ) one can show a number of properties. For instance, the concaveness of Cq^ 



for small m implies that for small m minimizing sequences are compact, and the minimizers 
are balls. One can also show that for sufficiently large m, a ball with volume m will be 
unstable with respect to symmetry-breaking perturbations; Ren and Wei have documented 
this phenomenon in two space dimensions |30] . 

For the three-dimensional case, however, one can show that balls become unstable with 
respect to splitting into two balls of half the volume before they become unstable with re- 
spect to small symmetry-breaking perturbations. This leads us to postulate the following 
characterization of global minimizers, when they exist: 



Conjecture 7.1. All global minimizers of the problem (4.1) are balls. 



Limiting structures. In both two and three dimensions, the limiting energies 'at the next 
level' Fq^ and Fq"^ penalize proximity of particles as if they were electrically charged. In two 



dimensions Lemma 6.2 guarantees that the masses m*, which play the role of the charges 
of the particles, are all the same; in three dimensions we conjecture that the same holds, 
although currently we are not able to exclude the possibility of (n — 1) equal masses and one 
different mass. 

The question whether minimizers of these Coulomb energies are necessarily periodic is 
a subtle one. It is easy to construct numerical examples of bounded domains on which 
minimizers can not be periodic; see e.g. [28] for examples on discs in R^. At the same time, 
the examples with many particles do show a tendency to a triangular packing away from the 
boundary. In the physical literature such structures are known as Wigner crystals, and in 
that field it is generally assumed that periodic structures have lowest energy. As far as we 
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know, there are no rigorous results that show periodicity without any a priori assumptions 
on the geometry. 

Turning to what can be proved, the closest related result we know is the two-dimensional. 



Leonard- Jones crystallization result of [Si]. Moreover, for the full problem (1.2), the only 
periodicity -like results we know of, in dimension larger than one, a statement concerning the 
uniformity of the energy distribution on large boxes [3] , and for finite-size structures in R" a 
scaling result bounding the energy in terms of lower- dimensional energies . 



The role of the mass constraint. Note that in the main theorems (6.1, |6.4[ 4.2, and 4.4) 



there is no mass constraint, as in (1.4), but only the weaker requirement that J v is bounded. 
This merits some remarks: 

• Free minimization of the limiting functionals Eq^, Fq^, Eg^, and Fq^ simply yields the 
zero function with zero energy. In order to have a non-trivial object in the limit some 
additional restriction is therefore necessary. Typically one expects to have a sequence Vrj 
for which the mass either is fixed or converges to a positive value. 

• The fact that only boundedness of J v is required also implies that this scaling of mass 
is the smallest one to give (for this scaling of the energies) non-trivial results; if / f 
converges to zero, then the limiting energies are also zero, and no structure can be 
determined. This conclusion resonates with the fact that in the formal phase diagram 



of the Ohta-Kawasaki functional (1.3) the phase at the extreme ends of the volume 



fraction range is the spherical phase [TD] . 



Related work. Our results are consistent with and complementary to two other recent 
studies in the regime of small volume fraction. In [29] Ren and Wei prove the existence 
of sphere-like solutions to the Euler-Lagrange equation of (1.1), and further investigate their 
stability. They also show that the centers of sphere-like solutions are close to global minimizers 
of an effective energy defined over delta measures which includes both a local energy defined 
over each point measure, and a Green's function interaction term which sets their location. 
While their results are similar in spirit to ours, they are based upon completely different 
techniques which are local rather than global. 

In |15[ Il9j the authors explore the dynamics of small spherical phases for a gradient flow 



for (1.2) with small volume fraction. Here one finds a separation of time scales for the dy- 
namics: Small particles both exchange material as in usual Ostwald ripening, and migrate 
because of an effectively repulsive nonlocal energetic term. Coarsening via mass diffusion 
only occurs while particle radii are small, and they eventually approach a finite equilibrium 
size. Migration, on the other hand, is responsible for producing self-organized patterns. By 
constructing approximations based upon an Ansatz of spherical particles similar to the clas- 
sical LSW (Lifshitz-Slyozov- Wagner) theory, one derives a finite dimensional dynamics for 
particle positions and radii. For large systems, kinetic- type equations which describe the evo- 
lution of a probability density are constructed. A separation of time scales between particle 
growth and migration allows for a variational characterization of spatially inhomogeneous 
quasi-equilibrium states. Heuristically this matches our findings of (a) a first order energy 
which is local and essentially driven by perimeter reduction, and (b) a Coulomb-like interac- 
tion energy, at the next level, responsible for placement and self organization of the pattern. It 
would be interesting if one could make these statements precise via the calculation of gradient 
flows and their connection with F-convergence |31j . 
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We further note that this asymptotic study has much in common with the asymptotic 
analysis of the well-known Ginzburg-Landau functional for the study of magnetic vortices (c/. 
|32| [TBI [2]). However our problem is much more direct as it pertains to the asymptotics of 
the support of minimizers. This is in strong contrast to the Ginzburg-Landau functional 
wherein one is concerned with an intrinsic vorticity quantity which is captured via a certain 
gauge-invariant Jacobian determinant of the order parameter. 



Although the energy functional (1.3) provides a relatively simple, and mathematically ac- 
cessible, description of patterns in this block copolymer system, rigorous results characterizing 
minimal-energy patterns in higher dimensions are few and far between. Apart from the work 
in this paper we should mention the uniform energy distribution results of [3] which provide 
a weak statement of uniformity in space, and the comparison of large, localized structures in 
multiple dimensions with extended lower-dimensional structures [H] . 

For a slightly a different model for block copolymer behavior additional results are avail- 
able. In |22^ l23l IM] the authors study an energy functional consisting of two terms as in £ 



in (1.2), but with the i?~^-norm replaced by the W~^'^-norm, or equivalently by the Wasser- 
stein distance of order 1. For this functional the authors study the symmetric regime, in 
which A and B appear in equal amounts; a parameter e characterizes the small length scale 
of the patterns. They prove that low-energy structures in two dimensions become increasingly 
stripe-like as e — > 0, that the stripe width approaches e, and that the Gamma-limit of the 
rescaled energy measures the square of the local stripe curvature. 
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